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Abstract 

Shear viscosity of the classical theory is measured using classical 
microcanonical simulation. To calculate the Kubo formula, we measure 
the energy-momentum tensor correlation function, and apply the Green- 
Kubo relation. Being a classical theory, the results depend on the cutoff 
which should be chosen in the range of the temperature. Comparison with 
experimentally accessible systems is also performed. 


1 Introduction 

Transport coefficients, in particular shear viscosity, are very hardly accessible 
quantities in perturbative quantum field theory calculations. Transport is char¬ 
acteristic to systems where information is spread by diffusion, the time evolution 
is ^ V Dt, the diffusion constant being the corresponding transport coefficient. 
The diffusion constant itself is proportional to the quasiparticle lifetime D ^ t. 
This is infinite in a free gas, and inversely proportional to some powers of the 
coupling constant for weak couplings. Therefore the perturbative evaluation of 
the Kubo formula ^ requires resummation of an infinite set of diagrams [2] • To 
circumvent this difficulty one can use effective methods to calculate the trans¬ 
port coefficients. One of these methods is the use of Boltzmann equations which 
is equivalent with the resummation of the singular part of the full perturbation 
series Hi]- Boltzmann equation method is used to obtain general results in 
gauge theories proving that to the leading order one has a shear viscos¬ 

ity r] ~ yring' Boltzmann equation methods are used also in other models to 
compute shear viscosity, like in meson models [71 [HI ini or in full QCD [TUI [11] • 
Other perturbation theory motivated methods to calculate the shear viscosity 
are 2PI resummation techniques [12] or the generalized quasiparticle approach 

m- 

Apart from the technical difficulties, also the applicability of perturbation 
theory makes these results less relevant for strongly interacting QCD-like sys¬ 
tems. Small value of the shear viscosity of the QCD plasma, reported by analy¬ 
ses of experimental data [14] suggests that the QCD matter is close to a perfect 
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liquid |15j . This implies that the interaction is rather strong, the quasiparticle 
lifetime is very short, and so perturbation theory is hardly applicable. 

Where perturbation theory is not well applicable, one seeks non-perturbative 
methods. Computer Monte Carlo (MC) simulation of QCD was used to extract 
shear viscosity data roughly in agreement with measurements m- The tem¬ 
poral range of the Euclidean formalism of the MC setup, however, makes the 
correlations less sensitive to long range physics which are relevant for trans¬ 
port m- Other popular method is to use the dual theory approach, based on 
AdS/CFT correspondence. Then weakly coupled five dimensional gravity can 
be used to compute transport coefficients in strongly coupled (conformal) field 
theories nails]. There are several model studies in this field which calculate 
shear viscosity by this method. 

Another nonperturbative method to approach the dynamics of quantum field 
theory is the use of classical theories to study both equilibrium [SO] [2T1 [22l 
Ha HH 123 and nonequilibrium phenomena [Sa HZl IMl UHl ■ Here one applies 
classical equations of motion starting from some initial conditions, and solve 
them by numerical methods on a finite mesh. The system thermalizesQ which 
in a classical system means equipartition of the energy. From the classical 
trajectories we can evaluate expectation values of different observables as time 
averages. 

From the point of view of perturbation theory classical and quantum systems 
are similar . In particular one can study expectation values of composite op¬ 
erators like (?/)). Comparing the classical and quantum computations 

one finds that with an appropriate choice of the cutoff of the classical theory 
Ac; T the quantum results can be nicely reproduced m- 

Encouraged by these results we tried to use classical simulations to compute 
the shear viscosity in a simple bosonic classical system, the model. The Kubo 
formula for the shear viscosity [T] contains commutator of spatial components of 
the energy-momentum tensor. We computed it with help of the Green-Kubo re¬ 
lation which is the classical counterpart of the quantum Kubo-Martin-Schwinger 
relation (fluctuation dissipation theorem) [JT]. This system has the potential 
to show a phase transition, similarly to the QCD case (although it is a second 
order here, as opposed to the crossover nature in QCD). This makes possible to 
study the rj/s ratio near the phase transition. 

The paper is organized as follows. First we overview the details of the dis¬ 
cretization and classical simulation method for the <I>^ model in Section II. In 
Section III we discuss the thermalization process and the measured character¬ 
istics of the thermal equilibrium, in particular thermal mass. In Section IV we 
report on our results of the energy-momentum tensor correlation functions and 
the classical values of the shear viscosity. In Section V we apply our method 
to quantum systems and present the 77 /s ratio, also in comparison with the 
experimentally measured values in different systems. The paper is closed with 
Conclusions. 

^Note that we work with finite systems with finite energy density where thermalization is 
possible. 


2 



2 The system, discretization and simulation al¬ 
gorithm 

The system we study is the quartic scalar model, which has the Hamiltonian 
density 

Hx = + i(V$(x))2 + + ^‘i>"(x). (1) 

Here $ denotes the held, H its canonical conjugate. The corresponding equations 
of motion (EoM) are 


$ = n, n = (2) 

6 

We remark that by rescaling the helds $/-\/A and H —>• n/-\/A, the equa¬ 
tions of motion become A-independent. We could work therefore with A = 1, 
but for better readability we keep the notation of A. 

We discretize the model on a symmetric hnite spacelike mesh 

3 

[7 = {x = ^ riiaei | = 0 ... — 1}, 

i=l 

where are orthogonal unit vectors and a is the lattice spacing; we express all 
dimensional quantities in lattice units and so we choose a = 1. For the lattice 
size we have in our simulations N = 36, 40 and 50, and use periodic boundary 
conditions. The discretized Laplacian is 

3 

A$(x) = ^ [$(x + e^) - 2$(x) H- $(x - e*)]. 

The discretized Hamiltonian can be written as H = X^xgc/ where the Hamil¬ 
tonian density formally equivalent to 0, "^ith (V$(x))2 = [$(x -h e,) - 4>(x)]^ . 

This is, however, not a local expression anymore, as it connects nearest neighbor 
held values. 

For the evaluation of expectation values we also need Fourier transformation. 

It is dehned on the reciprocal lattice U with the following dehnition (which 
corresponds to the fftw++ conventions |32j l 

/kGC7 = E /x, Uu = E h- (3) 

xGa kGd 


The reciprocal lattice is equivalent with the original lattice in case of cubic 
lattices we used. The Fourier-transformed Hamiltonian reads 


1 V- 

E 


kG(7 




24iV6 


E^ki4>k.$k3$k4- 


(4) 


kiG(7 
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where + X]i=i 4sin^ J2i ki = 0 in the last term. 

The time evolution in computer is realized using the leap-frog algorithm. 
Here one chooses a time step dt, so at nth step one arrives at time t = ndt. In 
the time step from n — 1 to n, the two equations of (H) are treated subsequently: 
first one evolves the held conhguration 

(5) 

then the canonically conjugated held conhguration, using the new values of the 
held: 

(x) = (x) +dt (x) - (x) - ^ (x)) , (6) 


with the discretized Laplacian. 

We can use the notion of the energy in the discretized model, too, as E = 
H = X^xGC/^' This quantity is conserved only for continuous time evolution; 
since we evolve the time in discrete steps, the total energy is not necessarily 
conserved. An important consistency check for the reliability of the algorithm 
is that in a long run the energy remains conserved. The leap-frog algorithm 
satishes this requirement. 

The classical ground state of the system, ie. the minimum of the energy is at 
spatially homogeneous held. If and A is positive, the minimum is reached at 
$ = n = 0. The coupling A must be positive otherwise the Hamiltonian of the 
system is not bounded from below. If < 0, the minimal energy is reached at 
a hnite |$| = <I>o value: this is the spontaneous symmetry broken (SSB) phase. 
The minimum condition yields <I>g = . 


3 Description of the thermal equilibrium 

Since we are primarily interested in the equilibrium properties of the system, we 
may start from an arbitrary initial condition. Practically we started with d>(x) = 
0, and from random values of n(x). After a certain time evolution (practically 
around t/a ^ 10000) we arrive at a steady equilibrium state. While the complete 
system forms a microcanonical ensemble, for local observables we can use a 
canonical ensemble. To determine its properties, we have taken the histogram 
of Hx, and found that it can be described by a Gaussian. This corresponds to 
the Boltzmann distribution ~ with /3 as a parameter interpreted as the 

“inverse temperature”. Therefore one can compute the expectation value of a 
local operator A(<i), H) which depends on the fields as 

(A($,n)) = ^ / n (7) 

d xeu 

where Z = J In Fourier space we should handle the 

problem that there is a relation between the integration variables 
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and similarly for Ilk, because of d) and 11 in coordinate space are real. To 
overcome this problem we introduce a purely real field 


~ _ J Re$k if fca > 0 
\ Im$k if ks < 0. 

This allows to write 

(A($,n)) = i f [] dnkdl>kA($,n)e-^^(‘^’n). 

ked 

To measure the temperature we use the relation 

OO 

(lEkP) = -^ / dflk (flk + n?.k) = 2N^T 

— OO 


( 8 ) 

(9) 


( 10 ) 


with T = 1//?. Using this formula we can check that the system arrived at 
equilibrium, by verifying that (lllkp) is independent of k (equipartition). An 
example for this distribution in a completely thermalized state is shown in Fig.[2 



Mode (A:) 


Figure 1: Mode dependence of temperature (fc^ = (27r^^)) 

We remark that the above canonical equilibrium description is in fact a 3D 
field theory of the initial conditions. As compared to the original action which 
is four dimensional, we have a dimensionally reduced theory. As a consequence 
the mass (energy) dimension of the field is [$] = 1/2. This fact will be used 
later when we apply dimensional analysis. 

In the thermal equilibrium we can perform perturbation theory. Although 
at large coupling (where we actually performed our simulations) results of per¬ 
turbation theory are not necessarily perfect, but in several aspects these may 
“guide the eye” to understand some robust features of the results. The details 
of perturbation theory in the classical theory can be found in [^. One uses 
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here two types of propagators 


Gret{k) = - 2 ’ iGsoik) = ^ g{k) ( 11 ) 

(fco + ley - ui% fco 

where the free spectral function is 

£>(fc) = 27rsgn(A:o)(5(fcg - w^), (12) 

and 

= ^4sin^ e 0 ... - 1. (13) 

As a final, technical issue, we remark that solving the field equation corre¬ 
sponds to the pure microcanonical, energy conserving approach to the thermo¬ 
dynamics. However, knowing that the system reaches equilibrium with Boltz¬ 
mann distribution, we can also use a canonical approach with a heat bath: in 
the language of the equations of motion it can be realized as a Langevin equa¬ 
tion. There we introduce a 7 damping parameter and a noise represented by 
^(x) independent stochastic variables with uniform distribution at each time 
step. We then change the update of H to 

n(”)(x) = (l- 7 dt)n("-i)(x)-hdt (^A$(”)(x) -to2$(”)(x) - ^ ($(”)(x)) Vc(x) 

(14) 

This stochastic process drives the system towards an equilibrium distribution 
with V{E) ^ distribution function. Because of the Einstein relation 

27 T = (^^) we can control the temperature of the thermal distribution. This 
algorithm can largely speed up the thermalization. After the system arrived at 
the equilibrium, we switched off the noise and damping terms in order that it 
does not influence the measurements. 


4 Equilibrium observables 

After we reached the equilibrium state, we can measure expectation values using 
time average 

to+t 

(A(ci>,n)) = i j (15) 

*0 

The equilibrium system can be characterized by a single value, for example the 
temperature. 

4.1 Energy 

We measured the relation between the temperature and the energy density, the 
results are shown in Fig. [^. We found that the relation is linear, with slightly 
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Figure 2: Temperature dependence of total energy {N^ = V = 50^), A = 5, 
= —0.5 

different slope in the symmetric and SSB regimes. One can clearly identify the 
phase transition region. Since our goal was not to study the phase transition 
point very accurately, we did not try to focus on this regime close enough to be 
able to tell details about it. 

This figure tells us that the heat capacity is proportional to the number of 
modes, as it is expected from a classical theory. Rewriting the lattice spacing 
a, this also means that the specific heat is proportional to a^. This is the well 
known Rayleigh instability (ultraviolet catastrophe) of the classical plasma. To 
have physically meaningful result, the lattice spacing must have a finite value. 

4.2 Mass and symmetry breaking 

It is important to notice, that the mass parameter of the Lagrangian (the bare 
mass) is not the same as the mass appearing in the observables (the effective 
mass). Physically it happens because of the nontrivial effect of the fluctuations. 

To estimate this effect (cf. Refs. [23 [57]) we used background field method. 
We shifted the classical field with its expectation value: $ —>■ $o + where 
{(fi) = 0. The shifted Lagrangian reads 

C = -—#2 _ Ai4 _ ( 2 ^ 1 2 _ ^2) _ A^2 2 _ 

2 24 4 ^ 

6 ^ 24^ ^ ^ 

To lowest order (Hartree approximation) we substitute the fluctuations by their 
expectation values. Using the fact that (ip) = 0 we find up to a constant 

47 ) 

This means that the effective mass is modified by the effect of the fluctuations. 
Since the mass dimension of the field is [p] = 1/2, by dimensional reasons 
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^ T. On the other hand this is a correction to the mass squared, and so 
the coefficient is also dimensionful, with finite lattice spacing it is proportional 
to a~^. The coefficient in leading order in perturbation theory reads 




N3 

kGU 


N—^cci 


Eli4sin2(5^)+m2 


777 T 

0.2527T-—+0(771^), (18) 
47r 


but this number is unreliable for large couplings. 

One consequence of this formula is that at fixed negative bare mass the 
effective mass term will be positive at high enough temperature. This means 
that the minimum of the effective action for the constant field (the constrained 
free energy) will become zero at this temperature: the symmetry is restored. 
To see it we measured the expectation value of the field for various tree level 
masses at different temperatures. The result can be seen on Fig. 
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Figure 3: Temperature dependence of the expectation value of the field. 

A more delicate question is that in the 3D classical field theory, unlike in 
the 4-dimensional theory, the temperature influences the renormalization. This 
means that the meaning of mass and temperature cannot be separated from 
each other in that clear way as it can be done in the 4D case. We therefore 
also performed simulations with fixed effective mass at different temperatures: 
this requires to tune the bare mass parameter. In practice we fix the desired 
effective mass and the bare mass, and tune the temperature accordingly with 
the application of Langevin equations described earlier. 

For the definition of the mass we measured the correlation function 

G(t,x):=($(t,x)$(0)) (19) 

In leading order of perturbation theory we expect that this correlator is the free 
one, where we should also take into account the mass modification: 

G(t,k) = (|4>fc|^) cosujkt, 
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( 20 ) 




where + m^ + ^ ( 93 ^). If one goes beyond the hrst order of perturbation 

theory, then one obtains self-energy corrections, and the pure harmonic behavior 
of the correlator will be spoiled. If a one-particle particle mass shell is dominant, 
then we can speak about quasiparticle excitations. In that case in real time 
evolution one can observe a damped oscillation 

G{t, k) ~ exp(—</rfc) cos(wfcf). (21) 

However, the closer we are to the phase transition point, the worse behavior 
could be observed for the static $ field, as it is demonstrated in Fig. We 




Figure 4: The real time behavior of the field correlation function at different 
bare masses, A = 5. 

can see that far from the phase transition point, where the effective mass is 
large, the quasiparticle assumption is valid. With decreasing effective mass the 
fit works worse and worse. In this case the definition of the notion “mass” is 
not unique anymore. For a more sophisticated description we should use the 
complete spectral function; but we just need a characterization of the mass and 
temperature. For that purpose we use the best quasiparticle fit to the real time 
data for the zero mode. This is some mean value of the spectral peak, in the 
vicinity of the phase transition point it remains finite, as opposed to the inverse 
spatial correlation length. 

The temperature dependence of the so-defined mass is shown in Fig. [^with 
fixed bare mass = —0.5. We can clearly see the position of the phase 
transition point which sits at the minimum of this curve. One may also check. 
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Figure 5: Temperature dependence of the effective mass parameter for various 
interaction strength (A) where = —0.5. It has a minimum at the phase 
transition point, but with our definition the minimum is not at zero. 


whether we reached the infinite volume (thermodynamical) limit. For that we 
determined the temperature dependence of the mass at various volumes, see left 
panel of Fig. This plot suggests that we reached already the thermodynamical 
limit. We can also check the temperature dependence of the effective mass, this 
is shown in the right panel of Fig. We see that for different bare masses the 
effective mass values sit on a unique curve. 
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Figure 6: Left panel: effective mass at different volumes. Right panel: effective 
mass in theories with different bare mass. 
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5 Viscosity 


The central topic of this paper is the determination of the shear viscosity. The 
Kubo-formula for momentum transport [I] requires to compute 


ri = lim 

w-fO 


QT12T21 tc 
UJ 


0 ) 


( 22 ) 


where 

0 AB(x) = ([^(a:),i?(O)]), (23) 

77 is the shear-viscosity and T 12 is the 12 component of the energy momentum 
tensor; in case of the scalar held theory it is T 12 = di^d 2 ^- In classical theory 
we cannot measure the commutator of two operators, but we can measure the 
correlation function instead. For A{x) and B{0) operators it is dehned as 


Sab{x) = {A{x)B{0))ci, 


(24) 


where the “cl” subscript refers to the classical correlation function. To connect 
this quantity with the viscosity we use the Green-Kubo formula, which claims 
that 

£'AB.c/(w, k) = /3ujSab,ci{uj, k), (25) 

which is a direct consequence of the quantum Kubo-Martin-Schwinger relation 

QABioJ, k) = (1 - e~^‘^)SAB,ci{uJ, k), (26) 

in the /3a; —> 0 limit. Using this relation the viscosity is 

Vci=PST,2T,2ik = 0). (27) 

The direct result of our simulations in real time can be seen in Fig. We 
repeated the measurements of S'ti 2 Ti 2 for hve different conhgurations, meaning 
that after each measurement we allowed the system evolve in time until we 
reached an independent configuration. This makes it possible to estimate the 
statistical error of the simulation. 

The relevant information, the transport-peak can be extracted from the 
Fourier-transformed data shown in Fig. To understand what we see in this 
figure we recall that in the leading order of perturbation theory we expect a 
branch cut starting at 2 to, and a Dirac-delta peak at k = 0, just like in the 
Fourier transform of (<i)^(a;)$^(0)). The higher order terms result in the smear¬ 
ing of the cut and the Dirac-delta peak as well, the former yielding a broad 
bump, the latter leading to the transport peak. The desired result is the height 
of this peak at a; = 0. 

One can repeat this measurement at different temperatures as well, this can 
be seen on Fig.|^ This figure suggests that the classical shear viscosity depends 
more or less linearly on the temperature. It also follows from dimensional anal¬ 
ysis: since [4>] = 1/2 is the mass dimension of the field, [T 12 ] ^ 3, its correlator 
has 6th power of energy dimension. After Fourier transformation there remains 
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Figure 7: T 12 autocorrelation function in real time. 
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Figure 8: T 12 correlator in frequency space. 
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2 , and after division by the temperature, there remains 1, a linear energy de¬ 
pendence. Since the main source of energy dependence in the classical case is 
the temperature, we expect proportionality with T. 

To verify numerically this, we also present the rjd/T curve in Fig. |10[ We 
see that the ratio is approximately constant, but with an enhanced behavior 
near the critical point. This figure suggests that the classical viscosity, just like 
other susceptibilities, shows a critical behavior, exhibiting a peak at the phase 
transition point. 

Classical viscosity 


0.8 
, 0.7 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Temperature (T) 


Figure 9: Temperature dependence of the classical shear-viscosity. 



Temperature (T) 


Figure 10: Classical shear viscosity over temperature as a function of the temper¬ 
ature. The enhanced behavior is around the phase transition point, suggesting 
that the viscosity also shows a critical behavior. 
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5.1 Interpretation 

One can compute the classical viscosity in perturbation theory to leading order, 
similarly what was done in Ref. |30] . With point splitting we can write 


S{x) = {Ti2{x)Ti2{0)) = hm d^id^'idyidy'i (l>(a:)l>(a:')l’(j/)l’(2/')) = 

y' -t-y—^O 

= d1iG3D{x)d2iG3D{x) + {did2iG3D{x)y. (28) 


After Fourier transformation we have 

Sip) =J ■^i;^[klip-k)l + kik2ip-k)i{p-k)2]iG3Dik)iG3D{p 

2 f d'^k klip - k)l + kik2ip - k)iip - k)2 n\ i i\ 

= ^ / oTAl- n -Td- Qik)gip-k). 

J koipo-ko) 


k) = 


(29) 


This is very similar to the correlation function discussed in 

very similar to the discontinuity of the quantum version of it. 


To proceed we use the Green-Kubo relation (271 to write 


f d^k kjkl 2n\ 

wv' ir® 


and also 


(30) 


Evaluating this expression with the free spectral function (12) we obtain infinity: 
this means that for free theories the viscosity, like all other transport coefficient, 
is infinite. We can apply a Breit-Wigner approximation for the spectral function 


Qik) 


4fco7fc 

{kl-u,)^+Akhl' 


(31) 


then we can approximate for small width and for k^ > 0: 


e^ik) 


27T 

IkUjl 


6iko - ujk). 


(32) 


This leads to 


Pel = 4r 


d^k klkl 

(27r)3 


(33) 


This integral is still divergent in the continuum limit. The leading contribution 
comes from large momenta. The asymptotics of 7 fcWfe has been found in |30j (see 
also [H]), it is . In the remaining factor one can put m = 0. Therefore 
to leading order we have 


15367r r klkl 

Pci ^2J' J ( 27 r )3 ( fc 2 _|_ ^(.2 _|_ ^ 2 ^ 3/2 


(34) 


Rewriting this integral in terms of the dimensionless momenta ka, and integrat¬ 
ing it in the Brillouin zone —tt < ka < tt, one finds 


_ Plat _ 321.37r 

Pci — ^4 ) Plat — ^2y 4" 


(35) 
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The fact that rjd ^ a~^ is consistent with the results of [30] for the <i>^ auto¬ 
correlation function. There a logarithmic divergence was found, therefore the 
energy-momentum tensor correlation function, which has two derivatives more 
than should scale as a~‘^. 

Although this result is just a first order perturbative estimate, the robust 
part of the result is that the classical shear viscosity is proportional to a“^. After 
determining rjiat by lattice simulations, this is the way how one can recover the 
classical value of the viscosity. 

The next question is how can one relate rjci to the shear viscosity of quantum 
systems. The facit of the comparison of the classical and quantum calculations 
[SOI [35] was that the perturbative quantum result was rather close to the classical 
one, if one chooses a cutoff A ^ T. The exact value of the coefficient is not known 
to be universal, probably it depends on the quantity in question. But, up to 
a constant, we can estimate the result of the full quantum result by setting 
r] ~ {aT)'^rid- If we measure the shear viscosity from the lattice, we have 

V ~ T^mat- (36) 


Below we will use unity for the proportionality constant. 

Of course, this formula is based on the assumption that the quantum result 
is dominated by the classical fields. Also, unfortunately, we do not know the 
coefficient in this formula, nevertheless, we can give a temperature profile of the 
77 /s ratio. 

The entropy density of the quantum system is poorly approximated by the 
classical modes, so we use for the ry/s estimate the entropy density of a free one 
component gas with 


s = 


uJk 

(27r)3 g/3<^fc _ X 


ln(l - 


(37) 


Here 


TTigjj, where rUe// was taken from the classical simulations. 


In this way we can present our estimate for the ry/s ratio in Fig. |1I[ We see a 
curve typical for the behavior of the shear viscosity in any matter near the phase 
transition point. For a qualitative comparison we show the ?y/s ratio for QCD 
[33] in Fig. 12 versus our results that was rescaled to fit to the high temperature 
part. 


6 Summary 

In this work we used numerical simulations of the quartic classical field theory 
to give an estimate for the shear viscosity. To this end we solved the discretized 
classical equations of motion. This leads to thermalization, where one can de¬ 
termine the expectation value of different observables by time averaging. We 
first determined the field autocorrelation function, and analyzed it by assuming 
quasiparticle behavior. Then we measured the correlator of the 12 component 
of the energy-momentum tensor. Using the Green-Kubo formula this quantity 


15 







4 

3.5 
3 

2.5 
CO 

er- 2 

1.5 
1 

0.5 
0 

0 0.2 0.4 0.6 0.8 1 1.2 

TjM 

Figure 11: The quantum rj/s estimated by the classical theory, using a = T~^ 
lattice spacing. 



% 


Lacey et al. 
l/4/7r 


-0.6 -0.4 -0.2 


Temperature (- 


Figure 12: Qualitative comparison of rj/s of QCD plasma and the classical 
theory, after rescaling. The QCD data are taken from |Mj. 
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is proportional to the shear viscosity. The shear viscosity ijiat was found to de¬ 
pend approximately linearly on T which was expected by classical dimensional 
analysis. This behaviour was superimposed by a characteristic critical behav¬ 
ior near the phase transition regime. Finally we pointed out that the classical 
viscosity comes from the lattice viscosity as rjci = 'rjiata~^ where a is the lattice 
spacing. Translating the classical result into the shear viscosity of the quantum 
system, we argued that ~ T is the correct choice, but the proportionality 
constant is not known. This allows us to make an estimate also on the temper¬ 
ature profile of rj/s, which turns out to be rather similar to the result of QCD 
near the critical region. 

As future prospects we plan to repeat this analysis to other models including 
gauge theories. Another interesting extension could be to study the effects of 
the quantum corrections to the equations of motion. This could give a hint on 
the reliability of the classical estimate. 


Acknowledgement 

The authors acknowledge useful discussions with A. Patkos, Zs. Szep. This 

research was supported by the grant K-104292 from the Hungarian Research 

Fund (OTKA). 

References 

[1] A. Hosoya, M. a. Sakagami and M. Takao, Annals Phys. 154, 229 (1984). 

[2] S. Jeon, Phys. Rev. D 52, 3591 (1995) [hep-ph/9409250]. 

[3] S. Jeon and L. G. Yaffe, Phys. Rev. D 53, 5799 (1996) [hep-ph/9512263]. 

[4] A. Jakovac, Phys. Rev. D 65, 085029 (2002) [hep-ph/0112188]. 

[5] P. B. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0011, 001 (2000) [hep- 
ph/0010177]. 

[6] P. B. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0305, 051 (2003) [hep- 
ph/0302165]. 

[7] A. Dobado and F. J. Llanes-Estrada, Phys. Rev. D 69, 116004 (2004) [hep- 
ph/0309324]. 

[8] M. Buballa, K. Heckmann and J. Wambach, Prog. Part. Nucl. Phys. 67, 348 
(2012) [arXiv: 1202.0724 [hep-ph]]. 

[9] S. Mitra and S. Sarkar, Phys. Rev. D 87, no. 9, 094026 (2013) 
[arXiv: 1303.6408 [hep-ph]]. 

[10] Z. Xu and G. Greiner, Phys. Rev. C 71, 064901 (2005) [hep-ph/0406278]. 


17 



[11] Z. Xu and C. Greiner, Phys. Rev. C 76, 024911 (2007) [hep-ph/0703233]. 

[12] G. Aarts and J. M. Martinez Resco, JHEP 0402, 061 (2004) [hep- 
ph/0402192]. 

[13] A. Peshier and W. Gassing, Phys. Rev. Lett. 94, 172301 (2005) [hep- 
ph/0502138]. 

[14] S. S. Adler et at. [PHENIX Gollaboration], Phys. Rev. Lett. 91, 182301 
(2003) [nucl-ex/0305013]. 

[15] J. Liao and V. Koch, Phys. Rev. C 81, 014902 (2010) [arXiv:0909.3105 
[hep-ph]]. 

[16] H. B. Meyer, Phys. Rev. D 76, 101701 (2007) [arXiv:0704.1801 [hep-lat]]. 

[17] P. Petreczky, J. Phys. G 35, 044033 (2008) [arXiv:0710.5561 [nucl-th]]. 

[18] P. Kovtun, D. T. Son and A. O. Starinets, JHEP 0310, 064 (2003) [hep- 
th/0309213]. 

[19] P. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 94, 111601 
(2005) [hep-th/0405231]. 

[20] G. Aarts and J. Smit, Phys. Lett. B 393, 395 (1997) [hep-ph/9610415]. 

[21] B. J. Nauta and C. G. van Weert, Phys. Lett. B 444, 463 (1998) [hep- 
ph/9709401]. 

[22] W. Buchmuller and A. Jakovac, Phys. Lett. B 407, 39 (1997) [hep- 
ph/9705452]. 

[23] W. Buchmuller and A. Jakovac, Nucl. Phys. B 521, 219 (1998) [hep- 
th/9712093]. 

[24] G. Aarts and J. Berges, Phys. Rev. Lett. 88, 041603 (2002) [hep- 
ph/0107129]. 

[25] B. Holdom, J. Phys. A 39, 7485 (2006) [hep-th/0602236]. 

[26] S. Borsanyi, A. Patkos and D. Sexty, Phys. Rev. D 66, 025014 (2002) 
[hep-ph/0203133]. 

[27] S. Borsanyi, A. Patkos and D. Sexty, Phys. Rev. D 68, 063512 (2003) 
[hep-ph/0303147]. 

[28] D. Boyanovsky, G. Destri and H. J. de Vega, Phys. Rev. D 69, 045003 
(2004) [hep-ph/0306124]. 

[29] J. Berges, B. Schenke, S. Schlichting and R. Venugopalan, Nucl. Phys. A 
931, 348 (2014) [arXiv:1409.1638 [hep-ph]]. 


18 



[30] A. Jakovac, Phys. Lett. B 446, 203 (1999) [hep-ph/9808349]. 

[31] M. Le Bellac, Thermal Field Theory (Cambridge University Press, 1996) 

[32] www.fftw.org 

[33] E. k. Wang, U. W. Heinz and X. f. Zhang, Phys. Rev. D 53, 5978 (1996) 
[hep-ph/9509331]. 

[34] R. A. Lacey, N. N. Ajitanand, J. M. Alexander, P. Chung, W. G. Holzmann, 
M. Issah, A. Taranenko and P. Danielewicz et al., Phys. Rev. Lett. 98, 092301 
(2007) [nucl-ex/0609025]. 


19 



